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Local Lorentz invariance violation can be realized by introducing extra tensor fields in the action 
that couple to matter. If the Lorentz violation is rotationally invariant in some frame, then it is 
characterized by an "aether", i.e. a unit timelike vector field. General covariance requires that the 
aether field be dynamical. In this paper we study the linearized theory of such an aether coupled to 
gravity and find the speeds and polarizations of all the wave modes in terms of the four constants 
appearing in the most general action at second order in derivatives. We find that in addition to the 
usual two transverse traceless metric modes, there are three coupled aether-metric modes. 



I. INTRODUCTION 

Recently there has been an explosion of research on the 
possibihty that Lorentz invariance is violated by quan- 
tum gravity effects (see e.g. ,1] and references therein). 
In a non-gravitational setting, it suffices to specify fixed 
background fields violating Lorentz symmetry in order to 
formulate the Lorentz violating (LV) matter dynamics. 
Hovifever, fixed background fields break general covari- 
ance. If we are to preserve the successes of general rela- 
tivity in accounting for gravitational phenomena, break- 
ing general covariance is not an option. The obvious 
alternative is to promote the LV background fields to dy- 
namical fields, governed by a generally covariant action. 
Virtually any configuration of any matter field breaks 
Lorentz invariance, but this differs in an important way 
from what we have in mind. The LV background fields 
we are contemplating are constrained either dynamically 
or kinematically not to vanish, so that every relevant 
field configuration violates local Lorentz symmetry ev- 
erywhere, even in the the "vacuum" . 

If the Lorentz violation preserves a three-dimensional 
rotation subgroup, then the background field must be 
only a timelike vector, which might be described by the 
gradient of a scalar, or by a vector field. In this paper we 
consider just the case where the LV field is a unit time- 
like vector which can be viewed as the minimal struc- 
ture required to determine a local preferred rest frame. 
We call this field the "aether", as it is ubiquitous and 
determines a locally preferred frame at every point of 
spacetime. Kinetic terms in the action couple the aether 
directly to the spacetime metric, in addition to any cou- 
plings that might be present between the aether and the 
matter fields. We refer to the system of the metric cou- 
pled to the aether as "Einstein-aether theory" . 

Here we investigate the linearized wave spectrum of 
this theory, and determine the complete set of mode 
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speeds and polarizations for generic values of the free 
parameters in the action. (Results for different special 
cases were previously published in Refs. These 
results identify the choices of constants in the action for 
which the linearized field equations are hyperbolic (and 
hence admit an initial value formulation), and they will 
be useful in extracting the observable consequences of 
such an aether field. 

Related work goes back at least to the 1970's, when 
Nordtvedt and Will began a study of vector-tensor theo- 
ries of gravity |4, 5 , 6, 7], which differed from the present 
work primarily in the fact that the norm of the vec- 
tor was not constrained. Gasperini, using a tetrad for- 
malism, studied in a series of papers 8] an equivalent 
formulation of the Einstein-aether theory studied here. 
Further related work has been done by Kostelecky and 
Samuel @ and Jacobson and Mattingly in the spe- 
cial case where the aether dynamics is Maxwell-like. The 
spherically symmetric weak field solutions were found for 
the gen eral Einstein-aether theory by Eling and Jacob- 
son [13. Vector-tensor theories have been studied in a 
cosmological context by Clayton and Moffatt ^H, J^J and 
Bassett et al The issues of causality and shocks in 
vector-tensor theories were studied by Clayton p!3. Fur- 
ther discussion on previous work can be found in 0, llOj . 
A proposal for Lorentz symmetry breaking via a scalar 
field with unusual kinetic term that makes the gradient 
tend to a timelike vector of constant norm has recently 
been investigated by Arkani-Hamed et al [Tsl flGf . Most 
recently, the issue of Lorentz violation in a gravitational 
setting has been examined in a systematic way by Kost- 
elecky [13. 



II. EINSTEIN-AETHER THEORY 

In the spirit of effective field theory, we consider 
a derivative expansion of the action for the metric 
gab and aether u°'. The most general action that is 
diffeomorphism-invariant and quadratic in derivatives is 

S = j d^X^(-R + Cu- Ku^'Ua - 1)) (1) 
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where 

Cu = -i^°^„„Va^/"Vfcu" (2) 

with 

R is the Ricci scalar, and A is a Lagrange multiplier 
that enforces the unit constraint. The metric signature 

is (H ), units are chosen such that c=l, and other 

than the signature choice we use the conventions of [l8l |. 

The presence of the Lagrange multiplier and C4 terms 
differentiate this theory from the vector-tensor theories 
considered in 0. The possible term RabU°'u'' is propor- 
tional to the difference of the C2 and C3 terms via integra- 
tion by parts and hence has been omitted. We have also 
omitted any matter coupling since we are interested here 
in the dynamics of the metric-aether sector in vacuum. 

Varying the action |^ with respect to w", g"'', and I 
yields the field equations 

Va - C4-UaV„i-u'' = \u„i (4) 

Gab = Tab (5) 
gabU^'u'' = 1 (6) 

where to compactify the notation we have defined 

rm = i^'^^™V^,u" (7) 

and 

u =u VqU , (8) 
and the aether stress tensor is 

Tab — ^ m{J(a'^Ui,') ~ J"^ (a^b) " J{ab)U"^) 

+Cl [(VmUaKV'^Ub) ~ iVaU„,)iVbU'^)] 
+ C4 ilaUb 

+ [Mn(V,„ J"*") - C4M^] UaUh 

-^QabCu- (9) 

In the above expression the constraint has been used to 
eliminate the term that arises from varying ^—g in the 
constraint term in J^l, and in the fourth line A has been 
eliminated using the aether field equation. 

A. Linearized field equations 

The first step in finding the wave modes is to linearize 
the field equations about the flat background solution 
with Minkowski metric rjab and constant unit vector u". 
The fields are expanded as 

gab = rjab + lab (10) 

u" = yt + v''. (11) 



The Lagrange multipler I vanishes in the background, 
so we use the same notation for the its linearized ver- 
sion. Indices are raised and lowered with rjab- We adopt 
Minkowski coordinates {x^^x^) aligned with u", i.e. for 
which rjab — diag(l, — 1, — 1, — 1) and — (1,0,0,0). 
The letters i, j, fc, I are reserved for spatial coordinate in- 
dices, and repeated spatial indices are summed with the 
Kronecker delta. 

Keeping only first order terms in u° and 7^6, the field 
equations become 

SaJ^'^^. = (12) 

Gil^ = T^a^ (13) 
«° + i7oo - (14) 

where the superscript (1) denotes the first order part of 
the corresponding quantity. The linearized Einstein ten- 
sor is 

41^ = -\^lab~\rab 

+7m(a,fc)" + \riab{Ul - 7mn ™"), (15) 

where 7 = 7m™ is the trace, while the linearized aether 
stress tensor is 

+K(a™j(i)™")]^i,u, (16) 

If we impose the linearized aether field equation 1)12(1 then 

the second and last terms of this expression for can- 
cel, yielding 

T^i.'^--5o./(%+9m^f:^'X), (17) 
The linearized quantity J^^^ is given by 

j'iib ^ ClV aUb+C2rjab'^mU"'- +C3VbUa+CiUa'^0Ub, (18) 

where the covariant derivatives of u°' are expanded to 
linear order, i.e. replaced by 

(Vatib)^^) = {Vb + 570fc),a + 57a6,0 - ^laO.b- (19) 

This completes an explicit display of the linearized field 
equations. 

The aether perturbations are coupled to metric per- 
turbations, due to the presence of the background aether 
vector u". Were it not for the aether background, the lin- 
earized aether stress tensor (|16|l would vanish, and the 
metric would drop out of the aether field equation, leav- 
ing all modes uncoupled. 

B. Gauge choice 

Diffeomorphism invariance of the action |^ implies 
that the field equations are tensorial, hence covariant 
under diffeomorphisms. The linearized equations inherit 
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the linearized version of this symmetry. To find the in- 
dependent physical wave modes we must fix the corre- 
sponding gauge symmetry. 

An infinitesimal diffeomorphism generated by a vector 
field transforms gab and m° by 

Sgab = ^i9ab = Va6 + Vfc^a, (20) 
Su" = r^M" = CV.^U" - u'^VmC- (21) 

In the linearized context, the vector field S,"' is itself first 
order in the perturbations, hence the linearized gauge 
transformations take the form 



lab 



= lab + da£.b + db(,a 



(22) 
(23) 



respect to time. Having made these choices of and ^i, 
holds, and the divergence of implies that 1(23 
holds. 

In the gauge (|24H25|I the tensors in the aether (|12ll and 
spatial metric equations (|13|l take the forms 



Jai,°' — Ci4(ui^00 — ^700,io) 



-C-lVikk 



Cl3lik,k0 



i% (07- 700,00 - IkLkl) 



(31) 



(32) 



Ti}^ = -ci3(w(ij)o + \li]..m) - kc2S^jjkk,oo (33) 
where we use the notation C14 := ci + C4, etc. 



The usual choice of gauge in vacuum GR is Lorentz 
gauge d''%b = 0, where 70b = jab - ^IVab- This gauge 
is chosen because it simplifies the Einstein tensor. The 
residual gauge freedom, which exploits the field equa- 
tions, further allows one to impose 70; — and 7 = 0. 
In the present case, the aether stress tensor ()17|l contains 
multiple terms in the derivatives of the metric pertur- 
bation and so Lorentz gauge is not particularly helpful. 
Moreover, the residual gauge freedom cannot be used to 
set 7oi and 7 to zero since these do not satisfy the wave 
equation. 

Instead, a convenient choice is to directly impose the 
four gauge conditions^ 



loi = 
ViA = 0. 



(24) 
(25) 



To see that this gauge is accessible, note that the gauge 
variations of 70^ and Vi^i are, according to (|22H23(I . 



S'-foi = Clq + ?o,i 



(26) 
(27) 



Thus to achieve the gauge H24l - I25|l we must choose ^0 and 
to satisfy equations of the form 



^1,0 + ^0,1 



Y. 



(28) 
(29) 



Subtracting the second equation from the divergence of 
the first gives 



^0, 



X,, 



(30) 



which determines ^0 up to constants of integration by 
solving Poisson's equation. Then can be determined 
up to a time-independent field by integrating H28|) with 



^ Alternatively, instead of setting Vi^i to zero it is equally conve- 
nient for finding the plane wave modes to set the spatial trace 
7ii to zero. 



III. WAVE MODES 

In General Relativity there are just two modes per spa- 
tial wave vector. Since u° has three independent degrees 
of freedom, we expect that in the Einstein-aether case 
there will be five modes all together. We now determine 
the wave modes in the chosen gauge. 

We assume a perturbation of the form 



lab 



^abe " 

^a^ikax" 



(34) 
(35) 



and choose coordinates such that the wavevector is 
(/co, 0, 0, /C3). The gauge conditions (|24H25|) then imply 



€3 



(36) 
(37) 



The problem is now to find the set of polarizations 
(cab, £a) and corresponding wave vectors ka for which the 
perturbation is a solution to the field equations H12l - I14|l . 

The component of the aether field equation H12I) is 
solved by definition of /, while the constraint equation 
(|14|l implies the relation 



(38) 



This leaves the spatial components of the aether equa- 
tion, together with the linearized Einstein equation. It 
suffices to use the spatial components of the Einstein 
equation, as the other components yield redundant in- 
formation (although they do provide useful algebraic 
checks). 

Inserting the plane wave ansatz ()34I35|) into the field 
equations yields 



[A3] 
[En] 



(ci4S^ 

ci4eoo 
£00 + 

[(1 + ' 



ci)ei 



2Ci3se/3 = 

Cl23e33 + C2e// = 

,2, I 1 
2 







C2)S 633 

ci23)s' - 1] eji = 



[En - E22] [(1 - ci3)s' - 1] (en - £22) = 



(39) 
(40) 

(41) 
(42) 
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TABLE I: Wave mode speeds and polarizations in the gauge 701 = Uj.i — 0. 



Wave mode squared speed s small c; limit polarization 

transverse, traceless metric 1/(1 — C13) 1 712, 711 = —722 

transverse aether (ci — ^cl + |c|)/ci4(l — C13) — > ci/cu 7/3 = [ci3/s(l — Ci3)]vi 

trace (ci23/ci4)(2 - C14)/ [2(1 + 02)^ - ci23(l + C2 + C123)] ^ C123/C14 700 = — 2uo 

711 = 722 = — Cl4l'0 

733 = (2C14/Ci23)(l + C2)vo 



[E12] [(1 - ci3)s2 - 1] = (43) 
[E13] ci3e/ + (-l + ci3)se/3 = (44) 
[£^33] (1 + C2)e// + £123633 = 0, (45) 

where [Ai] and [Eij] indicate the components of the 
aether and Einstein equations. We use the notation 
s = ko/k^ for the wave speed (which will be a true 
"speed" only when > 0), and the index / is dedi- 
cated to the two transverse spatial directions / = 1, 2, so 
that £// = en + £22 is the trace of the transverse spatial 
part of the metric polarization eab- 

We analyze the independent mode solutions assuming 
generic values of the constants 01^2,3,4- There are a total 
of five modes, two with an unexcited aether which corre- 
spond to the usual GR modes, two "transverse" aether- 
metric modes, and a fifth trace aether-metric mode. The 
two modes corresponding to the usual gravitational waves 
in GR are found when all polarization components van- 
ish except en, 622 and ei2. To avoid over-determining 
the speed s the trace equation [Ejj] must be identically 
satisfied, hence e// = 0. Then the [En — E22] and [E12] 
equations yield the speed. 

The two transverse aether-metric modes have nonzero 
polarization components e/ and £73, and the [Aj] and 
[£'73] equations together yield the speed and the ratio 
673/67. The fifth and final mode involves only 69 and 
the diagonal polarization components 6aa (no sum on 
a — 0,1, 2, 3). To avoid over-determining the speed, the 
difference equation [En — E22] — must be identically 
satisfied, hence 611 = 622 = ^ii/'^- Equations [A3] and 
[£'33] and the constraint (|38|) then allow all polarization 
components to be expressed in terms of just one, after 
which [-E77] determines the speed. The resulting mode 
polarizations and speeds are displayed in Table 

In the limit — > the transverse traceless modes 
become the usual gravitational waves, with unit speed. 
Note that these modes are entirely decoupled from the 
aether perturbations even when Ci ^ 0. 

The small Ci limits of the transverse aether and trace 
mode speeds depend on the ratios of the constants. If 
C2,3,4 vanish both speeds approach unity, but any other 
value is possible. The wave speeds and non-zero polar- 
ization components for the special case 02,3,4 = were 
previously reported in Q (the speed for the trace mode 
is inverted there), and the Maxwell-like case C13 = C2 = 
C4 = was analyzed in 2] (in both cases using different 



gauges). In the latter case, the transverse waves all have 
unit speed, while the trace mode has zero speed, so does 
not exist as a propagating wave. 

A peculiar special case occurs if C14 = 0, since then the 
aether wave speeds are generally infinite. This happens 
because no time derivatives of the aether field then arise 
in the field equation 131|) = 0. (The more special case 
Ci4 = C23 — was shown by Barbero and Villasehor [l9l | 
to be equivalent to general relativity via a m° dependent 
field redefinition of the metric.) 

When the constants Ci are chosen so that is positive 
and finite for all modes, the linearized equations are evi- 
dently hyperbolic. (It is not known whether this property 
extends to the nonlinear equations.) In these cases, since 
the dispersion relation = sfc is linear, |s| represents 
the signal propagation speed of disturbances. It is easily 
checked from Eqns. (41)-(45) that the Einstein tensor has 
nonvanishing components for each of the modes. (These 
equations display components of Gab — Tab, so just those 
of Gab remain when the Ci are set to zero.) Hence the 
modes all have gauge-invariant, physical significance. 

If is negative for a mode then the corresponding 
frequency is imaginary, indicating the existence of expo- 
nentially growing and decaying solutions. In such a case 
the theory is unstable and hence presumably unphysical. 



IV. OBSERVATIONAL APPLICATIONS 

An important open question is the sign of the energy of 
the various wave modes. To answer this, it is necessary to 
first determine the expression for energy in the linearized 
Einstein-aether theory, which has not yet been done. 

To compare the wave behavior of the theory with ob- 
servations the wave emission from astrophysical sources 
must be determined. To begin with, the analog of the 
quadrupole formula would enable the decay of binary 
pulsar orbits to be computed. Note that the presence 
of transverse aether and trace modes strongly suggests 
that dipole and monopole radiation will also exist and 
contribute to the energy loss. 

The wave emission depends of course on how the aether 
field couples to matter. A direct coupling could lead to 
local Lorentz violating effects which may exist but are al- 
ready quite constrained. However, even a small coupling 
to the matter source might be large enough to produce 
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an observable effect. Even without any direct coupling 
to matter, the extra modes will still be excited through 
their coupling to the time dependent metric produced by 
the moving matter sources. 

The results obtained here for the linearized theory 
should also be useful in computing the PPN parameters. 
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